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ABSTRACT 

Aims. We analyse multi-wavelength observations of 32 young star clusters and associations in M33 with known oxygen abundance 
(8 < 12 + log(0/H) < 8.7). The data set includes ultraviolet (UV), optical, mid-infrared (MIR), CO (1-0) and 21-cm line (Hi) 
observations. We derive the spectral energy distribution (SED) of these systems and the properties of their gaseous environment 
to investigate the process of star formation and the interplay with the interstellar medium (ISM). 

Methods. We determine age, bolometric luminosities, masses, and the extinction by comparing the multi-band integrated photometry 
to single-age stellar population models. The best-fit solutions have been obtained using the Large Magellanic Cloud (LMC) extinction 
curves. 

Results. The stellar system ages range between 2 and 15 Myr, masses are between 3 x 10 2 and 4x 10 4 M Q , and the intrinsic extinction, 
A v , varies from 0.3 to 1 mag. We find a correlation between age and extinction (young clusters being more reddened than older ones), 
and between the cluster mass and size. From the MIR emission we infer the presence of a dust component around the clusters whose 
fractional luminosity at 24 //m, L24/L So /, decreases with the galactocentric distance. However, the total infrared luminosity inferred 
from L 2 4 is smaller than what we derive from the extinction corrections. We also find that the Ho- luminosity predicted by population 
synthesis models is higher than the observed one, especially for low-mass systems (M < 10 4 M Q ). This difference is reduced, but not 
i-C ■ erased, when the incomplete sampling of the initial mass function (IMF) at the high-mass end is taken into account. 

Conclusions. Our results suggest that a non-negligible fraction of UV ionising and non-ionising radiation is leaking into the ISM 
outside the Hn regions. This agrees with the large UV and Ho- diffuse fractions observed in M33, but it implies that stellar systems 
younger than 3 Myr retain, on average, only 30% of their Lyman-continuum photons. However, the uncertainties in cluster ages and 
the stochastic fluctuations of the IMF do not allow us to accurately quantify this issue. We also consider the possibility that this 
discrepancy is the consequence of a suppressed or delayed formation of the most massive stars, given that it mainly affects the young 
and less massive clusters. We do not find any clear correlation between the cluster mass and the gas surface density or metallicity. 
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OO 1 . Introduction A number of studi es have analysed the cluster population of 

(N ■ M33 bo th from space dChandar et alj|1999t |2001|) and from the 

H !, m . . . • , ,. , . . . ground dMa et alJl200Ul2002ll2004l) 7aithough thev have mainly 

, • M33 is the nearest late-type spiral, its properties being in be- ° : , — : — ' — „ ■ ■ , 

V«0 1 * i- i vi »u »i-n ii r j in. j i concentrated on rather evolved stellar systems with ages between 

jST tween earlier-type spirals like the Milky Way and M3 1 and the ~ „ , . , . „ „ , , 

„ . , _c , • ■ , t iV^ , i i a few tens of Myr and tew Gyrs. Nonetheless, the population of 

O star-forming dwarf galaxies in the Local Group. With a relatively „ J , lri \. , . , r , . 

^.fil^L „™,-t a™ n a vm-3 vr. w-i t™.-2% y° un2 stellar s y stems (£ 10 Myr) in M33 has not been thor- 



high star-formation rate per unit area (3.4 x 10 M Q yr kpc ), , ? . . , \~ , / ,. , , , ,. , 
... \ v ■ J l inno k j i • v oughly investigated, and only a tew studies have been dedicated 

seven times higher than M31 (Kenmcutt 1998) and a low inch- ° i ^2 ■ 

to such young star-forming systems. 
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nation, M33 is an ideal environment to study the properties of 
young stellar systems and star-forming regions. 



A recent multi-wavelength analysis of radio-selected sources 
Young star clusters and associations are important for study- in M33 led to the identification of 1 1 young, optically-visible 
ing the ongoing star-formation activity and the chemical evo- stellar clusters withou t finding any highly obscured systems 
lution of their parent galaxy. The study of young stellar pop- (Buc kalew et al.l 120061) . Extinction seems generally lower in 
ulations enables us to analyse the initial phase of the cluster M33 than i n the Milky Way, and th e absence of large molecular 
evolution, which is mainly driven by hot massive stars. Despite complexes dRosolowskv et al]|2007l) may imply that even in an 
their short lifetime, these stars play an important role in galaxy early st ar-formation phase ste llar clusters may not be highly ex- 
evolution, because they are the main source of ionisation, heat- tincted. Corbel liet al.l ((2009) examined the properties of young 
ing, and the chemical enrichment of the interstellar medium stellar clusters in M33 selected at 24 pm with Ha counterparts 
(ISM). Massive stars trace the current star-formation activity, and stressed the importance of combining multiple wavelength 
they power the infrared (IR) radiation of dust and, at the end of observations (from the UV to MIR) to study star-forming sites 
their lives, they eventually trigger further star-formation events and to improve the und erstanding of massive star formation, 
via supernova explosions. Moreover, from the properties of the IVerlev et all (120071 12009) showed that MIR images of M33 ob- 
gas and dust distribution of the cluster environment one can bet- tained with the Spitzer telescope allow us to identify Hn regions 
ter understand the mutual relation between star formation and with 24 fim luminosities as low as 2 x 10 37 erg s~' and to study 
the ISM. the star-formation properties of the M33 disc. 
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Here we present a multi-wavelength analysis of a sample of 
32 young clusters and associations in M33 selected from th e cat- 
alogues of emission-line obj ects of iMagrini et al.l d2007bl) and 
Rosol owskv & Simonl d2008l) . Our aim was to define a list of 
young stellar systems selected on the basis of their Ha emis- 
sion, for which a measure of their metal abundance was avail- 
able. In addition, our objects have a 24 /mi counterpart in the 
catalogue of Verley et al. (2007). The final sample includes well 
known giant Hn regions such as NGC588 and IC132, located in 
the outskirts of the galaxy, large OB associations such as C400 
and B0222, and more compact Hn regions such as B0670 or 
LGCHII1 1 with sizes smaller than 10 pc. 

We use UV (GALEX), UBVRI (Local Group Survey), 
MIR (Spitzer IRAC and MIPS) observations to constrain ages, 
masses, extinction, and dust properties of the sample from the 
comparison to model spectral energy distributions (SEDs). We 
search for possible correlations among the derived parameters, 
we compare the properties of these stellar systems to their metal- 
licity, we look for radial gradients in the observed emissions at 
different wavelenghts, and finally, we use archival CO and Hi 
observations to relate the total stellar masses to the local gas 
content. 

The data are presented in Sect. 2, and in Sect. 3 we discuss 
the sample selection. The method to derive the photometric mea- 
surements in the different bands is described in Sect. 4, while an 
overview of the SED-fitting technique is given in Sect. 5. The 
results are discussed topic by topic in Sect. 6. The issue of the 
ionising photon deficiency within the clusters is analysed in Sect. 
7. The properties of the gaseous environment of the sample are 
presented in Sect. 8, and in Sect. 9 we summarise our results and 
conclusions. 

2. The data set: archival UV, optical and IR images 

2.1. Galex images 

We u sed Galaxy Evolution Explorer (GALEX; iMartin et al.l 
[2005b observations of M33 in the far-UV (FUV, A eff = 
1528 A) and near-UV (NUV, A eff = 2271 A), distributed by 
iGil de Paz etail (120071) . The angular resolution (FWHM) is ~ 
4" in the FUV and ~ 5".5 in the NUV images, corresponding to 
approximately 20 pc at the M33 distance. The exposure time was 
3400 sec in both filters. A more detailed description of GALEX 
observations of M33 and of the data reduction and calibration 
procedures can be found in Thilk er et al.l ((2005). 

2.2. Ground-based optical data: the Local Group Survey 

UBVRI images of M33 have been taken from the Local Group 
Surve y (LGS) of galaxies currently forming stars dMassev et al.l 
2006). M33 observations have been performed with the Kitt 
Peak National Observatory 4-m telescope using the wide-field 
Mosaic cameras. Narrowband images through the filters cen- 
tred on \Om] /1500 7, Ho- and [Sn] /1673 1 are also provided 
(Massev et al. 2007). The data are resampled to a scale of 0.27 
arcsec/pixel. The FWHM of the point spread function in the V 
band is l'.'l, corresponding to ~4 pc. More details about the ob- 
servations and data reduction can be found in these papers. 

We also made use of an optical narrowband Ha image of 
M33 obtained with the 0.6 meter Burrell-Schmidt telescope a t 
Kitt Peak National Observatory dHoopes & Walterbosl 12000). 
The 2048 x 2048 pixel CCD provides a field of view of 1.15 
x 1.15 square degrees, with a scale of 2.03 arcsec/pixel. The fi- 
nal mosaic covers an area of about 1.75 x 1.75 square degrees. 



Table 1. Sample of young clusters in M33 selected in this work. 
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32:43.0 
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25 


B0045 


01 


33:29.0 


30:40:24 


8.49 ± 0.04 " 


26 


B0029 


01 


33:47.8 


30:38:38 


8.21 ± 0.10* 


27 


BOO 16 


01 


33:50.1 


30:37:30 


7.98 ±0.19* 


28 


B0033b 


01 


33:34.9 


30:37:05 


8.40 ± 0.08 * 


29 


B0015b 


01 


33:50.3 


30:33:42 


8.73 ± 0.23 * 


30 


BOO 13c 


01 


33:54.1 


30:33:10 


8.26 ±0.12* 


31 


B0208f 


01 


33:34.1 


30:32:13 


8.07 ± 0.06 * 


32 
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01 


33:07.6 


30:31:01 
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Notes . The identi fi cation names are from: B-Boulest eix et al] dl974h : 
C - ICourtes et al.l (Ti987h : GDK9 94Gordon et alj d!999l) : LGC-HII- 
n4Magrini et al l J2007bl) : VGVC4Viallefond et alj 71989): MCM00 - 
Magri ni et alj d2000l) : MA 4Mavall & Allerl ( 119421) : S 4Searlel ( Ti97lh . 
Substructure designations within individual cl uster (indicated by th e let- 
ters at the end of their name) are taken from Hodge et al. 12002 |). Th e 
oxyg en abundance determinatio ns are from: (fl) Magri ni et al] d2007bh : 
^ iRosolowskv & Simonl J2008h 



More details on the observations and reduction procedure can be 
found in Hoopes & Walterbos (2000). 

2.3. Spitzer images 

Spitzer IRAC and MIPS images of M33 were obtained under 
the Guaranteed Time Observer (GTO) programme PID 5 (PI R. 
Gehrz). We used the IRAC (3.6-8 ^m) and MIPS (24 pm) maps 
of M33 produced bv lVerlev et al] d2007l) . The spatial resolution 
of the IRAC maps is between 1'.'7 and 2"(~ 7 pc), while at 24 pm 
it is 6"(24 pc). The map scale is 1.2 arcsec/pixel in both IRAC 
and MIPS images. 

2.4. CO and HI maps 

We used the 3 -mm CO map of M33 (J=l-0 rotational 
line) obtained from the combination of the Berkeley Illinois 
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Fig. 1. Location of the 32 clusters selected in this study, overlayed on the GALEX (FUV and NUV) images of M33. 



Maryland Association (BIMA) array (Engargi ola et alj [2003; 
Rosolows kv et alJ l2007) and the Fi ve College Radio A stronomy 
Observatory (FCRAO) single dish dHever et al.ll2004l) . The final 
resolution of the merged map is 13". The median noise level of 
the channel maps is 240 mK, corresponding to a mass sensitiv- 
ity of 5 xlO 4 M (4cr detection in two adjacent channels; see 
Rosolowsky et al. 2007) using a conversion factor X = 2 xlO 20 
cm" 2 (Kkms-'r 1 . 

The 21 -cm observations of M33 have been taken at 
the Westerbork Radio Syn thesis Telescope (WRST) array 
dDeul & van der HuTs11[l987l) . with a spatial resolution of 24" x 
48". 



3. Sample selection 

The sampl e was selected from t he catal ogues of emission line 
object s of Ma grini et al.l (l2007bl) and of iRosolowskv & Simonl 
(2008). The two catalogues contains 72 and 61 Hn regions re- 
spectively. The oxygen abundances were obtained using the 
temperature-sensitive emission line [Oin] /14363, with a direct 
measurement of the electron temperature. The condition of a 
measured and not assumed electron temperature is essential, be- 
cause the determination of the ionic and total chemical abun- 
dances exponentially dependends on it. 

The 32 Hn regions of the sample were selected on the basis 
of their shape and their position, preferentially round and/or lo- 
cated in regions of the galaxy not too crowded in the Ha map. 
Then we searched for a 24 jjm counterpart in the Spitzer/MlPS 
map of Verley et al. (2007), and by cross-correlating the two lists 
of objects we ended up with 25 sources with metal abundances 



within the range 8 < 12 + log(0/H) < 8.7. We then added seven 
more regions at large galactocentric radii, which have a clear de- 
tection in the 24 fim map even though they are not included in 
the Verley et al. (2007) catalogue. The final sample includes well 
known giant Hn regions such as NGC588 and IC132, which are 
located in the outskirts of the galaxy, large OB association as 
C400 and B0222, and more compact Hn regions such as B0670 
or LGCHII1 1 with sizes smaller than 10 pc. Figure Q] shows the 
position of the selected clusters and associations overlayed on 
the GALEX (FUV) image. In the rest of the paper we will refer 
to these young stellar systems by using the term "clusters" even 
though, in many cases, their size is larger than that of typical 
bound clusters. 



4. Multi-wavelength photometry 

We performed aperture photometry within circular regions 
choosing different apertures for each object. The aperture radius 
was defined on the basis of the extent of the V radial profiles, 
given the better resolution of the LGS images. The same posi- 
tion and aperture radius was applied to extract the photometry 
in the different bandpasses. The measurements were obtained 
using the IDL routine APER. Background subtraction was not 
straightforward because the clusters are found in different local 
environments. The contamination of nearby bright stars depends 
mainly on the object location in the galaxy and on the age of the 
background stellar population, and this is particulary critical in 
the R, I bands where the contamination by older stars is higher. 
The annuli for the background determination were selected indi- 
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Table 2. UV and optical photometric measurements. The last column displays the aperture radius used for the photometry. 



ID 


FUV 




NUV 




U 


B 


V 


R 




I 


' dp 




AB mag 


AB mag 


AB mag 


AB mag 


AB mag 


AB mag 


AB mag 


CO 


01 


15.22±0.00 


15. 


.31+0.00 


14 


,80+0.08 


14.63+0.06 


14.66+0.04 


14.76+0.04 


14 


.76+0.04 


11. 


02 


17.32±0.01 


17 


.40+0.01 


16 


,92+0.08 


16.95+0.06 


17.10+0.04 


17.18+0.04 


17 


.48+0.05 


4. 


03 


15.76±0.01 


15 


.69+0.00 


15 


,08+0.08 


14.87+0.06 


14.82+0.04 


15.24+0.04 


15 


.43+0.04 


9. 
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15.49±0.00 


15 


.29+0.00 


14 


,64+0.08 


14.71+0.06 


14.80+0.04 


14.92+0.04 


15 
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6. 
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14.84±0.00 


14 


.66+0.00 


14 


,59+0.08 


14.40+0.06 


14.44+0.04 


14.89+0.04 


15 


.16+0.04 


12. 


06 


15.20±0.01 


15 


.05+0.00 


14 


,86+0.08 


14.89+0.06 


14.95+0.03 


15.20+0.03 


15 


.26+0.04 


20. 


07 


14.42±0.01 


14 


.34+0.00 


14 


,09+0.08 


14.16+0.06 


14.30+0.03 


14.41+0.03 


14 


.64+0.06 


8. 


08 


16.23±0.04 


16 


.13+0.04 


15 


,67+0.08 


15.75+0.06 


15.90+0.04 


15.94+0.04 


16 


.22+0.05 


8. 
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18.67+0.11 
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.85+0.19 


17 


,61+0.08 


17.82+0.06 


17.97+0.04 


17.83+0.04 


18 


.34+0.06 


3. 


10 


19.28±0.04 


19 


.15+0.03 


18 


,36+0.08 


18.50+0.06 


18.52+0.04 


18.33+0.04 


18 


.60+0.06 


2. 


11 


18.69±0.03 


18 


.68+0.03 


17 


,51+0.08 


17.68+0.06 


17.80+0.05 


17.70+0.04 


17 


.94+0.06 


5. 


12 


14.85±0.00 


14 


.72+0.00 


14 


,38+0.08 


14.29+0.06 


14.32+0.04 


14.27+0.04 


13 


.97+0.04 


23. 


13 


16.74±0.01 


16 


.71+0.01 


16 


,13+0.08 


16.14+0.06 


16.24+0.04 


16.19+0.04 


16 


.20+0.04 


9. 


14 


17.02±0.01 


17 


.03+0.01 


16 


,66+0.08 


16.74+0.06 


16.87+0.04 


17.03+0.04 


17 


.34+0.05 


6. 


15 


18.07±0.02 


18 


.13+0.02 


17 


,20+0.08 


17.36+0.06 


17.50+0.04 


17.47+0.04 


17 


.81+0.05 


5. 


16 


19.67±0.09 


19 


.91+0.13 


18 


,82+0.08 


18.93+0.06 


19.01+0.05 


18.99+0.05 


19 


.30+0.08 


2. 


17 


17.29±0.02 


17 


.47+0.01 


17 


,16+0.08 


17.01+0.06 


17.05+0.04 


17.52+0.05 


17 


.65+0.06 


8. 


18 


18.39±0.03 


18 


.39+0.02 


18 


.19+0.08 


18.33+0.06 


18.46+0.04 


18.64+0.04 


18 


.96+0.05 


3. 


19 


19.28±0.04 


19 


.42+0.04 


17 


.90+0.08 


17.81+0.06 


17.76+0.04 


17.70+0.04 


17 


.94+0.05 


6. 


20 


19.90±0.04 


20 


.05+0.03 


19 


.46+0.08 


19.52+0.06 


19.55+0.04 


19.70+0.05 


20 


.19+0.08 


3. 


21 


16.89±0.01 


16 


.81+0.01 


16 


.60+0.08 


16.59+0.06 


16.72+0.04 


16.94+0.04 


17 


.17+0.05 


7. 


22 


18.37±0.02 


18 


.57+0.01 


17 


.67+0.08 


18.05+0.06 


18.24+0.04 


18.39+0.04 


18 


.59+0.07 


4. 


23 


18.09±0.03 


18 


.13+0.02 


17 


.27+0.08 


17.62+0.06 


17.77+0.04 


17.97+0.05 


18 


.15+0.08 


7. 


24 


19.33±0.04 


19 


,47+0.03 


18 


,48+0.08 


18.93+0.06 


19.11+0.05 


19.19+0.05 


19 


.29+0.09 


4. 


25 


16.25+0.01 


16 


.25+0.01 


15 


,58+0.08 


15.55+0.06 


15.66+0.04 


15. /2+0.04 


15 


,97+0.05 


9. 


26 


16.59±0.02 


16 


,60+0.02 


15 


,82+0.08 


15.91+0.06 


16.05+0.05 


16.07+0.05 


16 


,50+0.07 


7. 


27 


17.40±0.04 


17 


,44+0.01 


16 


,75+0.08 


16.94+0.07 


17.17+0.07 


17.22+0.06 


17 


,54+0.09 


5. 


28 


17.97±0.02 


17 


.89+0.02 


16 


,61+0.08 


16.68+0.06 


16.76+0.04 


16.62+0.04 


16 


,86+0.05 


6. 


29 


17.63±0.03 


17 


.65+0.03 


16 


.57+0.08 


16.75+0.07 


16.84+0.05 


16.69+0.05 


16 


,84+0.09 


9. 


30 


16.23±0.04 


16 


,02+0.03 


15 


,12+0.08 


15.21+0.06 


15.25+0.04 


15.09+0.04 


15 


,17+0.05 


11. 


31 


16.78±0.01 


16 


,63+0.01 


15 


,99+0.08 


16.09+0.06 


16.30+0.04 


16.36+0.04 


16 


,74+0.06 


7. 


32 


17.95±0.05 


17. 


,82+0.05 


16 


,94+0.08 


16.88+0.06 


17.01+0.05 


17.10+0.05 


17. 


,27+0.06 


7. 



vidually for each cluster in regions where the contamination of 
bright nearby objects was minimal. 

Table[2]displays the UV and optical photometry of the cluster 
sample and the corresponding uncertainties, which include the 
errors from the object count rate, the sky variance, and the zero 
point calibration. 

4.1. UV photometry 

The flux density in the FUV and NUV filters was measured 
within circular apertures given in Table [2l and the median back- 
ground flux was estimated in annuli of 5 pixel width, cho- 
sen in a way to reduce the contamination of nearby sources. 
Conversion from count rates to the AB photometric system 
(IOkeTll990) was obtained using the zero-points ZP^^ = 18.82 
and ZP^y = 20.08 (iMorrissev et al.ll2005l) . UV magnitudes 
were corrected for foreground Galactic extinction following 
Schlegel et al.l dl998h before starting the analysis of the cluster 
SEDs. 

4.2. UBVRI and Ha photometry 

UBVRI instrumental magnitudes were converted t o the standard 
Johnso n and Cousin system using stars from the iMassev et al.l 
(2006) catalogue. The colour term in the UVBRI transformation 
(described in Massey et al. 2006) was ignored, as it is smaller 
than the errors we report. The contribution of the Galactic fore- 



ground extinction was also subtracted from the optical magni- 
tudes according to Schleg el et al.l (11998). 

The continuum emission in the Ha image was removed by 
subtracting the /?-band image scaled by a factor 0.40, which was 
determined from the aperture photometry of 20 field stars in both 
filters. The calibration factor of 1.79 xl0~ 19 erg s _1 cirT 2 given 
in Massey et al. (2007) was applied to the count rates measured 
within apertures that included the total Ha emission around the 
cluster. The fluxes are displayed in Table [3] As a further check, 
Ha fluxes were determined for the same apertures on the map of 
Hoopes & Walterbos (2000), adopting a calibration factor of 8.2 
xlO -18 erg s crrT 2 . The resulting measurements agreed with 
the LGS map within 5% (lcr). 

Ha and oxygen emission lines may significantly contribute 
to the total flux in the R, V and U bands. We estimated the con- 
tamination from [Oin] /15007 and Ha lines to the photometry in 
the V and R filters, using the narrowband images provided by 
LGS. The correction was obtained by subtracting the flux mea- 
sured in the narrowband image within the same aperture used for 
the optical photometry (last column in Table 2). We used the cal- 
ibration factor of 2.1 1 xlO" 18 and 6.41 xl0~ I8 ergs~' cirT 2 A -1 
for the Ha and [Oin] images respectively, as given in Massey et 
al. (2007). The U photometry is affected by the contamination 
from the [On] A3121 line, but the LGS does not provide a nar- 
rowband image centred on this wavelength. Thus we multiplied 
the [Oin] narrowband fluxes for the TOnl /TOm l line r atio derived 
for the clusters in the iRosolowskv & Simon! (120081) catalogue. 
For those without a direct measurement of this ratio, we applied 
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Table 3. Ha and Spitzer IRAC/MIPS photometry of the sample. Flux units are in ergs s _1 cm"' 



ID 


r Ha 


XFf 




XFf 




AFf 




10- 13 cgs 


cgs 


10~ 12 cgs 


10~ 12 cgs 


10~ 12 cgs 


10~ 12 cgs 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 



25. 

4. 
52. 
26. 
61. 

4. 
32. 
11. 

4. 

2. 

2. 
19. 

7. 

3. 

3. 

0. 

5. 

0. 

4. 

0. 

4. 

0. 

1. 

0. 
25. 
15. 

4. 
12. 

6. 
16. 
18. 

1. 



29±0.13 
82±0.15 
85±1.62 
33±0.56 
11±0.36 
80±0.24 
64±0.38 
48±0.34 
18±0.09 
51 ±0.07 
58±0.06 
92±0.13 
29±0.05 
15±0.08 
66±0.07 
72±0.02 
94±0.09 
72±0.02 
13±0.08 
30±0.01 
88+0.15 
68+0.01 
78+0.04 
43±0.01 
56±0.87 
36±0.34 
14±0.14 
52±0.25 
74±0.23 
18±0.53 
67±0.47 
53+0.06 



6.10±0.74 
0.22±0.10 
1.20±0.68 
2.65±0.38 
1.57±0.26 
1.92±0.32 
1.92±0.31 
4.08±0.55 
0.42±0.12 
0.28±0.09 
0.32±0.10 
11.96±1.35 
1.62±0.28 
0.54±0.14 
0.42±0.12 
0.08+0.05 
0.12±0.06 
0.10±0.05 
0.16±0.07 
0.03±0.03 
0.17±0.08 
0.10±0.05 
0.09+0.05 
0.08+0.07 
4.76±0.60 
1.34±0.26 
1.54±0.26 
1.93±0.31 
0.91±0.21 
3.32±0.46 
1.78±0.29 
0.44+0.12 



4.23±0.53 
0.17±0.07 
1.00±0.37 
2.16+0.31 
1.37±0.22 
1.06±0.20 
1.42±0.23 
2.51+0.35 
0.30±0.09 
0.20±0.07 
0.22±0.07 
7.63±0.88 
1.20±0.21 
0.37±0.10 
0.37±0.10 
0.05±0.04 
0.10±0.05 
0.06±0.04 
0.12+0.05 
0.02±0.02 
0.12±0.06 
0.06±0.04 
0.07±0.04 
0.05±0.04 
3.96±0.50 
0.92±0.19 
1.09±0.19 
1.55±0.25 
0.63±0.15 
2.09±0.31 
1.11+0.20 
0.28±0.09 



9.73±1.06 
0.50±0.11 
1.54±0.27 
6.43±0.73 
2.31+0.31 
3.07±0.39 
3.90±0.47 

10.67+1.16 
1.26±0.20 
0.76±0.14 
0.93±0.16 

18.95+1.98 
3.12±0.39 
1.42±0.21 
1.18+0.19 
0.21±0.06 
0.16±0.06 
0.15±0.05 
0.29±0.08 
0.07±0.04 
0.18±0.06 
0.15±0.05 
0.08+0.04 
0.10±0.04 

14.34+1.52 
5.04±0.60 
5.15±0.60 
4.39±0.52 
3. 11 ±0.40 
9.32±1.03 
6.18±0.70 
1.42±0.21 



24.18+2.51 
1.65±0.23 
7.76±0.88 
16.40+1.71 
9.63±1.03 
7.04±0.77 
26.65±2.75 
26.15±2.69 
5.10±0.64 
3. 31 ±0.40 
2.76±0.35 
35.46±3.61 
7.98±0.87 
3. 81 ±0.44 
3.29±0.39 
1.21+0.17 
2.11±0.28 
0.46±0.09 
2.14±0.27 
0.23±0.06 
1.30+0.21 
0.40±0.08 
0.48±0.10 
0.23±0.06 
39.06±3.99 
34.17±3.60 
16.57+1.73 
12.29+1.33 
4.64±0.66 
28.48±2.96 
45.39±4.63 
3.29±0.39 



23.78±2.40 
0.75±0.09 
6.71±0.69 

15.77±1.60 
9.85+1.01 
3.14+0.33 

18.56+1.88 

16.17+1.64 
2.20+0.24 
2.39+0.26 
0.75+0.10 

24.89+2.51 
4.32+0.45 
1.89+0.21 
5.69+0.59 
0.38+0.06 
0.58+0.08 
0.19+0.03 
0.46+0.06 
0.09+0.02 
0.51+0.08 
0.15+0.03 
0.13+0.03 
0.07+0.02 

72.34+7.26 
8.75+0.91 
8.75+0.90 

12.10+1.23 
2.70+0.30 

11.66+1.19 

12.18+1.24 
2.20+0.24 



an average value of [On] /[Oiii] = 1 . In some cases this correction 
was not sufficient to reduce the excess in the U -band photometry. 

4.3. IRAC and MIPS (24 p.m) photometry 

Photometry at 3.6, 4.5, and 5.8 fim was derived through the same 
apertures used for the UV-optical bands, while the larger aper- 
tures chosen for the Ha photometry were adopted to measure 8 
and 24 fim fluxes. Uncertainties include Poisson noise photon 
flux, background fluctuati ons and a 10% flux calibration uncer- 
tainty for both th e IRAC dFazio et al.|[ 2004: Re ach et al.ll2005l) 
and MIPS bands dEngelbracht et al.ll2007l) ~Tablel3l displays the 
fluxes at 3.6, 4.5, 5.8, 8.0, and 24 ytzm for the sample. 



5. SED fitting with CHORIZOS 

CHi-square cOde for parameteRized modelling and characteri- 
zation of photometry and Spectrophotometry (CHORIZOS) is 
a x 2 minimization code that enables one to fit an arbitrary fam- 
ily of spectral energy di stributions to a set of photometric and 
spectrophotometric data (Mafz-Apellaniz 2004). It compares an 
observed set of magnitudes to a grid of model SEDs and cal- 
culates the probability for a given model to be consistent with 
the observed photometry. It allows one to constrain both in- 
trinsic parameters (such as age and metallicity) and extrinsic 
ones such as extinction law and colour excess. The user selects 
the SED family and chooses the behaviour of each parameter 
to be fixed, constrained within a certain range, or to be kept 



free. The code uses the monochromatic equivalents of E(B - V) 
and R v , £(4405-5495) and ^5495, where 4405 and 5495 are the 
assumed central wavelengths in Angstrom of the B and V fil- 
ter s. It allows on e to ch oose between the galactic extinction law 
oflCardelliet al.1 (1 19891) . the average LMC and LMC2 laws of 
Miss elt et all d 19991) and the Small Magel lanic Curve (SMC) ex- 
tinction law of [Gordon & Clavtonl (Il998l) . Then ax 2 analysis is 
performed to assign a likelihood to each SED as a function of 
the range of parameters. 

The code includes a set of STARBURST99 dLeifherer et al.l 
119991) instantaneous burst models for a ran ge of ages and metal- 
licities, produced using a Salpeter (1955) stellar initial mass 
function (IMF) with masses within the range 0.1-100 M . The 
metallicity of the cluster models was chosen at Z=0.004, 0.008 
and 0.02, while the age varied between 1 and 100 Myr. We ex- 
tinguished the theoretical SEDs with all the available extinction 
laws using a routine provided with the code. 

We ran CHORIZOS for all the objects in our sample using 
seven magnitudes (from FUV to / bands). We did not include 
the 3.6 and 4.5 /mi IRAC bands in the fitting becau se the con- 
tribution of hot dust emiss ion starts at A >3 /urn (jCarico et al. 
[T98ilKennicutt et alj2003l) and the STA RBURST99 models be - 
come inapplicable above that wavelength dLeifherer et al.ll999h . 
A more detailed analysis of the MIR region of the cluster SEDs 
will be presented elsewhere (Giovanardi et al. 2010, in prepara- 
tion). 

We performed the fitting procedure by setting the metallicity 
at the values displayed in Table 1, leaving the other parameters 
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X (A) X (A) X (A) 

Fig. 2. Output best-fit SEDs in units of AB magnitudes, overlayed on the measured integrated magnitudes from 
Vertical bars indicate the uncertainty on the photometry, while the horizontal ones show the approximate width 
numbers at the top-left corners of the plots correspond to the cluster ID given in TableQ] 



X (A) 

1000 to 10000 A. 
of each filter. The 



(i.e. age and intrinsic extinction in M33) unconstrained. Then 
we reran the code a second time, with a higher resolution grid, 
to refine the estimate of the parameters. The errors on the age 
and the extinction are estimated by the procedure. To assess the 
goodness of the fit, we relied on the value of x 1 provided by 
CHORIZOS; on average clusters with^- 2 < 10 were considered 
as good fits. 

Figures |2] and [3] show our SED-fitting results. For each clus- 
ter we display the output best-fit SED in units of AB magnitudes 
from 1000 to 10000 A overlaid on our measurements to show 
the quality of the fitting (Fig. [2]), and the model SED in flux units 
(A¥ A , erg s -1 cm 2 ) from 1000 to 3 xlO 5 A to compare the UV, 



optical and IR emission of the clusters (Fig. [3]). The numbers at 
the top-left and bottom-right corners of the plots correspond to 
the cluster ID given in Table [T] The best-fit SEDs overlaid on 
the observed magnitudes (Fig. |2]l show that in some cases the 
U photometry disagrees with the model, despite our attempt to 
take into account the effects of [On] line contamination, causing 
an increase in the value of^ 2 . Comparison between the best- fit 
model SEDs and the MIR photometry (Fig.O shows an excess 
emission at 3.6 fim and 4.5 fim bands, suggesting that this non- 
stellar emission may be due to a hot dust component. 
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Fig. 3. Model SED for each cluster in flux units (AF A , erg s _1 cm 2 ) from 1000 to 3 x 10 s A compared to the Spitzer IRAC and MIPS 
photometry. The numbers at the bottom-right corners of the plots correspond to the cluster ID given in TableQ] 



6. Results: comparison of cluster properties 

6.1. Ages, masses, sizes, and colours 

The SED fittings confirm the young age of the clusters, being 
younger than ~15 Myr. The results of the fitting procedure are 
displayed in Table [4] Uncertainties in age are given in Table [4] 
and are less than 0.2 (in log years). Variations in the IMF slope 
from the standard Salpeter value, or stochastic effects in the IMF 
when sampling small mass clusters, may increase the quoted un- 



certainties dCervino & Luridianal200l:ICorbelli et al.ll2009h . We 
will discuss the issue of stochastic fluctuations in Sect. 7.1. 

CHORIZOS also determines the bolometric correction to the 
V magnitude for each best-fit model. Knowing the age, metal- 
licity, and the bolometric luminosity, we inferred the stellar 
mass using the STARBURST99 evolutionary synthesis models. 
Figure[4](fo/7-Ze/f panel) displays the distribution in age and mass 
of the clusters. It shows that most of the clusters in our sample 
are found in a very narrow range of ages, while the mass distri- 
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Table 4. Results from the SED-fitting technique. 



ID 


log(age) 


A 


V 


E(B - V) 


log(L B „/) 


log(M) 




[yr] 


mag 


mag 


[L Q ] 


[M Q ] 


1 


6.97±0.01 


0.28:i 


:0.01 


0.13+0.01 


6.71+0.01 


4.31+0.01 


2 


6.51±0.14 


0.58:1 


:0.11 


0.26±0.03 


6.39+0.08 


3.24+0.08 


3 


6.68±0.09 


0.68:1 


:0.11 


0.30+0.05 


7.23+0.01 


4.37+0.01 


4 


6.72±0.03 


0.59:1 


:0.07 


0.24±0.03 


6.98+0.05 


4.20+0.05 


5 


6.54±0.03 


0.61=t 


:0.02 


0.23±0.01 


7.56+0.01 


4.49+0.01 


6 


6.65±0.15 


0.52:1 


:0.17 


0.23±0.07 


6.76+0.12 


3.86+0.12 


7 


6.71 ±0.02 


0.37:1 


:0.04 


0.14±0.02 


7.11+0.03 


4.31+0.03 


8 


6.48±0.16 


0.70:1 


:0.09 


0.29±0.03 


6.94+0.07 


3.78+0.07 


9 


6.40±0.10 


0.82d 


:0.07 


0.35±0.02 


6.21+0.03 


3.07+0.03 


10 


6.35±0.21 


0.93d 


:0.09 


0.39+0.02 


6.00+0.10 


2.88+0.10 


11 


6.36±0.13 


0.93d 


:0.06 


0.41±0.01 


6.37+0.04 


3.24+0.04 


12 


6.99±0.04 


0.40:1 


:0.02 


0.16+0.01 


6.95+0.01 


4.57+0.01 


13 


7.18±0.04 


0.39d 


:0.02 


0.16+0.01 


6.19+0.01 


4.04+0.01 


14 


6.44±0.16 


0.61:1 


:0.08 


0.25+0.02 


6.52+0.07 


3.37+0.07 


15 


6.38±0.16 


0.82:1 


:0.07 


0.36+0.02 


6.51+0.05 


3.38+0.05 


16 


6.36±0.15 


0.83d 


:0.07 


0.37+0.02 


5.85+0.04 


2.72+0.04 


17 


6.35±0.13 


0.56:1 


:0.05 


0.26+0.01 


6.44+0.03 


3.32+0.03 


18 


6.43±0.17 


0.54±0.08 


0.22+0.02 


5.89+0.07 


2.75+0.07 


19 


6.22±0.12 


1.05:1 


:0.03 


0.51+0.01 


6.49+0.01 


3.39+0.01 


20 


6.35±0.14 


0.66:1 


:0.06 


0.30+0.01 


5.57+0.04 


2.45+0.04 


21 


6.53±0.16 


0.56:1 


:0.10 


0.24+0.04 


6.46+0.09 


3.36+0.09 


22 


6.37±0.09 


0.61d 


:0.03 


0.28+0.01 


6.11+0.02 


2.98+0.02 


23 


6.43±0.13 


0.71d 


:0.08 


0.30+0.02 


6.28+0.06 


3.14+0.06 


24 


6.39±0.10 


0.66:1 


:0.05 


0.30+0.01 


5.75+0.03 


2.62+0.03 


Id 


£. A O i r\ 1 A 

6.48+0. 14 


0.70:1 


:0.12 


U.32+U.U3 


on i n (yi 
6.W+0.0 / 


1 O Q i n (Yl 

3.83+0.0 / 


26 


6.42±0.14 


0.80:1 


:0.07 


0.34+0.02 


7.01+0.06 


3.87+0.06 


27 


6.50±0.17 


0.66:1 


:0.11 


0.30+0.03 


6.50+0.09 


3.33+0.09 


28 


6.33±0.14 


1.05d 


:0.07 


0.46+0.01 


6.84+0.03 


3.72+0.03 


29 


6.20±0.11 


0.88:1 


:0.06 


0.40+0.01 


6.71+0.03 


3.58+0.03 


30 


7.16+0.1 1 


0.56:1 


:0.11 


0.27+0.09 


6.60+0.02 


4.44+0.02 


31 


6.46±0.17 


0.81:1 


:0.08 


0.33+0.03 


6.90+0.07 


3.75+0.07 


32 


6.49±0.14 


0.88:1 


:0.11 


0.38+0.03 


6.57+0.07 


3.40+0.07 



bution encompasses two orders of magnitude, between ~ 3 x 10 2 
and ~ 4xl0 4 M . In this and in the following figures we used 
different symbols (and colours) for clusters within different age 
bins: filled (blue) dots correspond to ages younger than 3 Myr, 
(green) squares to clusters with ages in the range 3-8 Myr, while 
the oldest ones (> 8 Myr) are indicated with (red) triangles. 

Observations of Giant Molecular Clouds (GMCs) in dif- 
ferent environments indicate that t hey follow a mass - radius 
relation M oc R 2 dLarsonl Il981h IWilson & Scovilld 119901; 
lHarris & Pudritzlll994l) . Thus, once we derived the masses, we 
checked whether the stellar systems under study follow a simi- 
lar mass-radius relation or not. As an estimate of the size of the 
cluster we used the half-light radius measured in the V band. As 
we show in the upper- right panel of Fig. [4] we find a correlation 
between the half-light radii of the clusters and the mass obtained 
from the SED fitting, with a slope of 1.5 + 0.2. The correlation 
is stronger for clusters with sizes smaller than 10 pc, while the 
others show a larger scatter around the best-fit line. 

Compact star clusters (R < 10 pc) with ages between few 
tens and few hun dreds of Myr show a poor corre lation between 
mass and radius dZepf et alj|1999t iLarsenl 12004). implying that 
the GMC mass-radius relatio n must be brok en at a certain stage 
during the cluster evolution. lLarse n (2004) find M oc R 0A in a 
sample of relatively young compact clusters (age < 1 Gyr and 
R < 10 pc) in nearby spiral galaxies. On the other hand, cluster 
complexes with a range of ages more similar to our sample (< 
10 Myr) but larger sizes (R > 85 pc ) follow a similar mass- 
radius relation to GMCs, M oc R 2 3 (Basti an et ai]l2005h . Our 



results seem to indicate that even at scales smaller than those 
of stellar complexes young stellar systems still display a similar 
trend to the GMCs, and that the observed stellar densities are 
representative of the proto-cluster cloud environment. 

The bottom-left panel of Fig. [4] shows the bolometric lumi- 
nosity distribution of the clusters, obtained from the SED fitting. 
The plot of the dereddened (B - V)o colour as a function of age 
(bottom-right panel) gives an overview of the colour range of the 
sample. Given their young age, the clusters are characterised by 
very blue colours and, as expected, "older" clusters are redder 
than the younger ones. 

6.2. Metallicity 

The metallicity gradient in M33, from 1 to ~7 k pc is rather shal- 
low dVilchez et alJll988b iMagrini et all l2007bl) . with a central 
abundance of 12 +log(Q/H) =8.36+ 0.04 and a slope of 0.027+ 
0.012 dex kpc~' (Rosolowsky & Simon 2008). This can be ex- 
plained by the global chemical evo lution of the M33 di sc through 
an inside-out formation process (Mag rini et alJl20O7ah . In Fig. 
[5] (top panel), the metallicity of the clusters is plotted versus 
the galactocentric distance, and it is shown that the clusters fol- 
low the radial trend found by Rosolowsky and Simon (2008) 
indicated by the solid line. However, Rosolowsky and Simon 
(2008) find a scatter of about 0. 1 dex in metallicity at a given ra- 
dius, which is unexplained b y the abundance uncertainties only 
dRosolowskv & Sim on 2008). A possible explanation would be 
that self-enrichment could be more efficient in more massive 
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Fig. 4. Top-left: cluster ages and masses determined from the SED fitting. Top-right: mass versus radius relation: the solid line with 
a slope of 1 .5 provides the best linear fits to the data points. The dotted-line shows the minimum size corresponding to the resolution 
of the V-band image. Bottom-left: bolometric luminosity distribution of the sample. Bottom-right: (B - V)q colour distribution as 
a function of age. The symbols correspond to three different age bins: log(age) < 6.5 (filled dots), 6.5 < log(age) < 6.9 (squares), 
log(age) > 6.9 (triangles). 



clusters, and produce the scatter around the average value. If that 
is the case, we should find a relation between mass and metallic - 
ity in our sample. 

The central panel of Fig. shows that this correlation is 
lacking in our sample, indicating that the metal abundance is re- 
lated to the local environmental properties of the clusters rather 
than being influenced by the cluster mass. Neither do we find 
a clear trend between the metallicity and age (lower panel); the 
Spearman correlation coefficient is -0.32, indicating a poor de- 
gree of correlation between these two parameters. 

6.3. Extinction 

Massey & Hutchings (1983) first concluded from the analy- 
sis of the IUE spectra of M33 Hn regions that its extinc- 
tion curve is different from the galactic one, and that it is 
characterised by a weak 2175 A hump similar to that of the 



LMC. iDodorico & Benvenutil (119831) find a similar result after 
a study of the Hn region IC132. For the brightest Hn region of 
M33, NGC604, a reddening param eter E(B-V)=0.12 has been 
obtain ed assuming a LMC curve (Gonzalez Delga do & Perezl 
I2000h . 

The E(B - V) values we derive for the M33 clusters range be- 
tween 0.1 and 0.5 mag, corresponding to a visual extinction Ay 
between 0.3 and 1 mag. Chandar et al. (1999) obtain E(B - V) 
values between 0.06 and 0.33 mag , with the majo rity of the clus- 
ters being found around 0. 1 mag. lPark & Led d2007l) derived red- 
dening values between 0.05 and 0.20 mag. Yet both samples are 
different from ours because they have a large range of ages as 
mentioned in Sect. 1. The LMC average and LMC2 curves of 
Misselt et al. (1999) were the ones providing the best fits in the 
majority of the cases, even though the SMC extinction curve was 
adopted for one cluster (C400). Additional evidence supporting 
the choice of either the LMC average or the LMC2 extinction 
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curve comes from the UV colours. The SMC extinction law pre- 
dicts larger variations in the FUV-NUV colour than what was 
observed in our sample. 

Figure |6]shows that there is a correlation between the cluster 
age and Ay (top-left). Younger clusters are more reddened, be- 
cause they are still partially embedded in the parent molecular 
cloud where they formed. Ay does not appear to clearly corre- 
late with the metallicity (top-centre), and neither does it correlate 
with the 8 pm (top-right) nor the 24 pm luminosities (bottom- 
left) normalised to the bolometric one. A trend between the ex- 
tinction and the total hydrogen column density is also lacking 
(H2+H1 +Hn ; see also Sect. 8), as displayed in the bottom-right 
panel. 

6.4. Radial gradients 

We searched for possible variations of the cluster properties with 
the galactocentric radius, especially with regard to their UV, Ha 
or MIR emissions (see Fig. [7j. 

Neither Ay (top-left) nor the L/7£/y/Lg / ratio (not displayed) 
show a clear trend with the galactocentric distance. A radial 
change in Ay would imply a decreasing dust at tenuation, as it 
is found for example in M5 1 ( Calzetti et al.l2 005). In this galaxy 
the UV colour, Ay, and the age of the clusters do vary with 
the distance, with the clusters becoming younger and less ex- 
tinct as one moves to the outer regions of the galaxy. Extinction 
is related to the dust-to-gas ratio, which depends on the metal- 
licity of a galaxy. Because steep metallicity gradients are ob- 
served in many galaxies, one would expect that the attenuation 
properties would also change with the gal actocentric distance . 
The shallow metallicity grad ient in M33 (iMagrini et al.l 12007a: 
IRosolowskv & Simonl 120 08) might explain the lack of a clear 
correlation between Ay, orLfw /Lb,,i, and the galactocentric ra- 
dius. 

The L° h JLboI ratio (top-centre, Ha corrected for extinction) 
decreases with galactocentric distance for the youngest clusters 
(filled dots). However, since the youngest age bin includes sev- 
eral low-mass clusters, this might be due to a selection effect: 



usters and associations in M33 

most of the selected low-mass clusters are at large galactocen- 
tric radii in regions that are not too crowded in the Ha map (see 
the last panel in the figure). We will discuss this in more detail 
in Sect. 7. Note that the Ha luminosity of the clusters with ages 
around 10 Myr is more than one order of magnitude higher than 
what STARBURST99 predicts, implying that subsequent star- 
formation episodes must have occurred in order to explain the 
observed Ha emission. 

We see a well defined radial trend for the MIR fractional 
luminosities. Both the I^/L^,,/ (top-right) and the Lg/L° FUV ra- 
tios (bottom-left) decrease by at least one order of magnitude 
at large galactocentric radii, at the edge of the star-forming 
disc. The MIR emission at 24 pm is mainly due to very small 
grain s with effective radii between 15 and 40 A dDraine & Lil 
l2007h . Variations in the abundance of small grains with respect 
to big ones at larger radii, or, more generally, a radial gradient 
in the dust-to gas ratio or in the dust tempera ture, may also con- 
tribu te to the decrease of the 24 nm emi ssion dHabartetalJ200U 
lAbereel et al.l2002tlCaizetti et al.l2005h . The luminosity at 8 pm 
is likely dominated by Polyciclic Aromatic Hydrocarbons (PAH) 
features around star-forming regions excited by the cluster UV 
radiation, and we find a good correlation between 8 pm and 
the FUV luminosities of the clusters (Fig. [8]). The PAH emis- 
sion is reduced in low-metallicity objects, but over the small 
range of metal abundances in our sample there is no clear depen- 
dence of the LsfL F uv ratio on the metallicity. At larger radii the 
volume density of the ISM decreases, therefore one would ex- 
pect that dust and PAHs may be more efficiently removed from 
star-f orming sites in the outer regions of the disc (Roussel et al. 
l200l . 

The total stellar mass does not seem to be clearly related to 
the cluster position within the galaxy as shown in the bottom- 
right panel of Fig. [7] even though the smallest mass clusters of 
the sample (M < 10 3 M Q ) are mainly located at large galacto- 
centric distances. 

6.5. Bolometric luminosity 

The SED-fitting method provides the value of the bolomet- 
ric luminosity for each cluster, Lgoi (see Table 5). We com- 
pared this parameter with the bolometric luminosities estimated 
from observations of the total IR and UV emission, i.e. L„ = 
Lfuv + Lnuv + L7vs/(1 - Co), where co ~ 0.5 is the albedo 
for a standard dust compo sition and grain size distribution 
dWeingartner & Draindl200ll) . 

The sum of the FUV and NUV luminosities provides a mea- 
sure of the cluster emission below 3000 A, which is close to the 
bolometric lumi nosity if the cluster is you ng and unobscured by 
dust dBell et al.ll2005l iThilker et alll2007h . If extinction correc- 
tions are unknown, one can use Ltir to estimate the fraction of 
the UV emission absorbed by the dust. Therefore Ut is consid- 
ered to provide a good appro ximation of the bolometric lumi- 
nositi es of young clusters (see lThilker et alJfeOOTblCorbelli et al.l 
120091) . Verley et al. (2009) have shown that the total IR luminos- 
ity in M33 can be estimated from the 24 pm emission and it is 
approximately Ltir ~ 12 x L24 (see Eq. 4 in Verley et al. 2009). 

Figure [9] shows that the model bolometric luminosities are 
higher than the observed ones for < 3 x 10 39 erg s^ 1 . This 
difference could be due to two main reasons: a) the assumption 
of a complete IMF may be incorrect when we consider clusters 
with bolometric luminosity (and mass) below a certain thresh- 
old, leading to an overestimate of the total cluster luminosity; b) 
the estimated total IR emission of the clusters is lower than the 
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amount of energy missing in the non-ionising UV region of the 
spectrum, implying that there is a deficit of non-ionising radia- 
tion in young clusters with bolometric luminosities lower than 
3 x 10 39 erg s . We will discuss these issues in more detail in 
the next section. 



7. The deficit of ionising radiation 

In the previous section we pointed out that the total IR emission 
of the clusters with log(LBo/) < 39.5 is lower than the amount 
of energy missing in the non-ionising UV region of the spec- 
trum. Here we want to compare the observed Ha luminosity with 
that predicted by an instantaneous burst model to check whether 
there is a deficit of UV-ionising radiation in the sample of clus- 
ters under study. 

We ran STARBUST99 to obtain the expected Ha luminosity 
of each cluster at the age, metallicity, and mass we determined, 
and we plot in Fig. QJjjthe ratio of the extinction-corrected ob- 
served Ha luminosity to the model one - L^ ff /L^ 99 - versus the 
age (left-panel) and the mass (central-panel) of the clusters. The 
figures show that this ratio decreases with age and mass, while 
the right-panel shows a shallow radial trend with the galacto- 
centric radius for the youngest clusters. There is a significant 
discrepancy between the observed and model Ha emission es- 
pecially in young (t < 3 Myr) and low-mass (log(M/M ) < 3.5) 



systems, for which L^/Ltj,* 99 is less than 0.5, or -0.3 in logarith- 



mic units. This would imply that more than 50% of the ionising 
flux is missing from these clusters, a higher percentage than what 
was found in other studies (see Sect. 7.3). 

The error bars in Fig.[10]show how the ratio of Ha luminosi- 
ties varies if one takes into account the uncertainties of the best- 
fit ages (Table 4). The upper errors bars are larger for clusters 
older than 2.5-3 Myr, because the model Ha emission rapidly 
decreases after this epoch. For the youngest clusters (t < 3 Myr), 
the mean L° Ha /L s H B J >9 ratio is 0.19*q^ 2 , while for those with an 
age between 3 and 8 Myr the mean is L° Ha /L s H B a " = 0.65+°* 3 . 
On the other hand, clusters older than 8 Myr (triangles in Fig. 
[TOt are very luminous at Ha wavelengths for their age, because 
L° Ha is more than 10 times higher than the prediction of simple 
stellar population models. It is likely that in these cases the as- 
sumption of an instantaneous burst is incorrect and an extended 
star-formation activity is required to explain the observed Ha 
emission. 

Hence, even taking into account the uncertainties of the clus- 
ter ages, the missing fraction of ionising radiation is rather large 
in young clusters, and we examine four different hypotheses to 
explain this issue. 

First we discuss how the incomplete sampling of the IMF 
and statistical fluctuations can affect the prediction of cluster lu- 
minosities. Then we consider the effects related to the cluster 
environment such as dust absorption and photon-leakage (i.e. a 
large fraction of UV radiation escapes the Hn region). Finally we 
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the same as for Fig. [4] 



briefly discuss whether the delayed formation of the most mas- 
sive stars or a truncated IMF can offer alternative explanations 
for the low Ha luminosity observed in very young clusters. 

7. 1 . The consequence of sampling effects in the IMF 

The library of models that we used for the SED fitting assumes a 
standard Salpeter IMF, fully populated up to a high-mass end of 
100 M . It is possible that the difference between the expected 
and observed Ha luminosities in the youngest clusters is because 
we are overestimating the number of massive stars owing to the 
incompleteness of the IMF at the high-mass end. 

Incomplete sampling is often neglected, and it is assumed 
that the mass spectrum of stellar clusters can be represented by a 
continuous function given by the underlying IMF. On the other 
hand, numerical simulations of synthetic clusters by a stochas- 
tic generation of stars show that this assumption is correct for 
massive systems (M > 10 4 M ), wherea s large deviations a t 
the high-mass end occur for small clusters dCorbelli et all 2009). 
The critical mass below which statistical fluctuations become 
relevant depends on the clus ter age and metallicity, a nd on the 
spectral range unde r study (|Cervino et alJl2 002. 2003]). For ex- 
ample, according to IFatuzzo & Ad ams (2008), sampling effects 
dominate in the FUV range when extinction corrected FUV 
fluxes are below 10 39 erg s _1 . Most of our selected objects lie 



above this threshold, but the situation might change if we are 
interested in the UV ionising radiation, because it has a steeper 
dependence on the stellar mass than the non-ionising continuum. 
Hence, to estimate the Ha luminosities of the smallest clusters 
in our sample, sampling effects need to be considered. 

In Fig.Q~T|we show how the median bolometric (diamonds) 
and Ha (stars) luminosities of a cluster d epend on its mass. Th e 
simulation was carried out as described in Corbe lli et alJ {2009), 
but here we assumed a maximum stellar mass of 100 M to be 
consistent with the fitting procedure used in this paper. The lu- 
minosities are computed for zero-age main-sequence stars. 

Low-mass clusters have a larger dispersion around the me- 
dian Ha luminosity than around the bolometric luminosity. 
Moreover, L//„ deviates considerably from a simple linear rela- 
tion (solid line) at masses lower than 1450 M . Below this limit 
the IMF is not fully populated at its high-end, and the stronger 
dependence of on the stellar mass implies that the mean 
value of Lboi/Lhc increases, despite the occasional presence of 
outliers. 

The IMF statistical fluctuations introduce a large disper- 
sion in the observable properties of the low-mass clusters (M 
< 1450 M ), when the linear relation between cluster mass and 
Lh„ breaks down. Thus, to account for the incomplete sampling 
of the IMF in these clusters, we will use the Ha luminosity and 
the dispersions derived with our simulation (see Fig. fTTb . Note 
that all clusters with M < 1450 M in our sample are quite 
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young, hence their expected Ha emission is well represented 
by the simulated data in Fig. QT| For older and massive clus- 
ters we shall use the luminosities obtained for a complete IMF 
(STARBURST99) w ith the dispersion given by our simulation. 
ICervino et al .] d20Q3h have shown that the minimum mass be- 
low which we expect stochastic effects, as well as the dispersion 
around the mean or median value, have a very weak dependence 
on the cluster age during the first 10 Myr. 

The expected Ha luminosities of low-mass clusters are lower 
when stochastic effects are considered, as shown in Fig. [12] The 
error bars in the figure include both the dispersion of the model 
Ha luminosity and the uncertainty due to the error on the age 
of the clusters. The mean L° Ha /L™° a d raticQ for the youngest clus- 
ters (f < 3Myr) is 0.30^q'q 7 . Therefore, taking into account the 



stochastic effects reduces the discrepancy between the expected 
and observed values, but the large uncertainties prevent us to 
draw any firm conclusions on how the missing Ha flux depends 
on the cluster mass, age, or galactocentric radius, as we show in 

Rg. EH 

7.2. Dust absorption 

Dust within or around the cluster can absorb the UV photons 
emitted by the hot massive stars, reducing the amount of ionising 
radiation, and contributing to the lack of observed Ha photons. 

The good spatial correlation between 24 /mi and Ha emis- 
sions shows that the dust is predominantly heated by OB stars 
within the Hn regions. The Spitzer resolution at 70 and 160 /mi 
is not sufficient to resolve our sources, and we estimated the to- 
tal IR luminosity from the 24 /mi emission following Verley et 
al. (2009). This agrees with what other studies found when the 
IR luminosity is dominated by emission at wavelengths longer 
than 24 /mi. Nevertheless, the total IR luminosities are lower 
than what is required by the modest extinctions obtained from 
the SED-fitting procedure. Some additional contribution to Ljir 



1 Note that to derive the mean ratio we discarded the two lowest mass 
clusters with the largest error bars. 
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Fig. 9. Comparison between the cluster bolometric luminosity 
determined from the SED spectral fitting (y-axis) and the ob- 
served bolometric luminosity derived by adding the measured 
FUV, NUV, and the total IR luminosity estimated following 
Verley et al. (2009) as Ltir ~ 12x L24. Symbols are the same as 
for Fig. H 



might come from a hot dust component heated by young stars in 
these sources. The SED indicates an excess above the expected 
stellar emission at 3.6 and 4.5 /mi. However, given the limited 
MIR emission detected around the clusters and, more generally, 
the low dust content of M33, it is difficult to explain the miss- 
ing ionising radiation as an effect of dust absorption only (see 
Sects. 6.4 and 6.5), even though this mechanism would explain 
the observed trend with cluster age. 

On the other hand, the way dust is distributed within or 
around the Hn regions may produce different attenuation of 
the emitted light, and the usual assumption of a homogeneous 
dust screen ma y lead to an underestimation of the effects of 
dust extinction dPetersen & G ammelgaard 1997). For example, 
a clumpy dust distribution would explain the missing fraction 
of ionising radiatio n (40%) observed in the super star clusters 
of SBS 0335-052 dReines et al.ll2008l) . A more detailed analy- 
sis of the dust properties for these clusters, taking into account 
more complex dust distributions, will be addressed elsewhere 
(Giovanardi et al. 2010, in preparation). 

7.3. Leakage of UV photons 

UV photons not absorbed by dust grains can escape from the 
star-forming region and explain the origin of the extended 
diffuse UV emissio n observed by GALEX in the M33 disc 
dThilker et al.ll2005h . 

The importance of taking into account the leakage of ion- 
ising photons from Hn regions has been increasingly dis- 
cussed in recent years, especially after the discovery of the 
diffuse ionised gas (DIG) in the ISM of the Milky Way 
(Reynolds et al.ll 19741). and in the discs of other spiral galaxies 
dGreenawalt et al.l ll998). The DIG components of galaxies such 
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Fig. 10. Ratio of the observed (extinction corrected) to the model Ha luminosity determined with STARBURST99 vs. the age of 
the clusters (left-panel), the mass of the clusters (central-panel), and the galactocentric distance for clusters with log(age) < 6.5 
(right-panel). The error bars give the range of the Ha luminosity predicted by STARBURST99 for each cluster due to uncertainty 
on the age. Symbols are the same as for Fig. [4] 



as M31 dWalterbos & Braunll 19961). t he LM C (lOev & Kennicuttl 
[19971 119981) . the SMC dRelano et alj|2002h . amount to the 30% 
and 5 0% of the total Ha luminosity. St udies of field OB stars in 
M33 dHoopes & Walterbosl2000l l2003l) have given evidence that 
at least 20% - 30% Lyman continuum photons must be able to 
escape from the Hn regions in this galaxy. Nonetheless, there is 
still a debate about whether the energy budget of photons leaking 
out of Hn regions is sufficient to io nise the DIG or if additional 
mechanisms should be considered dVoges et alJ l2008). 

Theoretical models have shown that the amount of leak- 
age is regulated by the structure and du mpiness of the ISM 
(Mill er & Coxll993l:IWood & Mafhisl2004l) . Inhomogeneities in 
the medium favour the leakage of radiation because the presence 
of low-d ensity regions may increase the mea n free path of the 
photons dVoges et alJ2008l:lReines et alJ2008l) . The leakage also 
depends on the spectrum of the ionising radiation, and on the 
evolutionary stage of the cluster, with the fraction of escaping 
photons being larger when stars experience a hard Wolf-Rayet 
phase, which occurs, depe nding on the mass and m etallicity of a 
star, between 3 and 6 Myr (Castel lanos et al.ll2002l) . 

Relating the leakage of ionising photons to the beginning of 
the W-R phase (Castellanos et al. 2002) may explain the age- 
dependence of the L^ ff /L^ 99 ratio for clusters older than 3 Myr, 
but the issue would still remain open for the youngest systems. 

7.4. Deficiency or delayed formation of massive stars 

The previous considerations cast some doubts on leakage and 
dust absorption as the only processes causing the deficiency of 
ionising radiation. 

If the star-formation process is not instantaneous, clusters 
might reach their expected Ha luminosity only after all stars 
have formed, especially the most massive ones, which determine 
the cluster Ha luminosity rather than the total mass. Massive 
stars might form later during the cluster evolution either because 
they n eed particu lar physical conditions and processes to take 
place dPallall2005l) or because the time-scale for the gas to col- 
lapse is longer in a low -density environment such as the disc 
of M33 dElmeg reen 1999). Therefore, the delayed formation of 
the most massive stars appears as an interesting and promising 
way to explain why only very young clusters are underluminous 
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Fig. 11. Expected median bolometric (diamonds) and Ha (stars) 
luminosity of stellar clusters as a function of cluster mass and 
their corresponding uncertainties. 

in Ha, and this should be investigated in detail in future work. 
Particularly, it would be interesting to know whether the excess 
of Ha emission in the oldest clusters could also be explained by 
this scenario without the assumption of multiple burst events. 

A strong deficiency of massive stars compared to what 
is predicted by a stochastically sampled IMF should also be 
considered. For example, truncated models of the stellar IMF 
(iMeurer et al.l 120091: iPflamm-Altenburg et al] 120091 and refer- 
ences therein) predict lower cluster luminosities than stochas- 
tically sampled models (Corbelli et al. 2009). This is because 
truncated models link the cluster mass to the most massive stars 
that can form in a cluster, without considering the occasional 
presence of outliers. The IMF over the whole galaxy in this case 
would deviate from the power law assumed in each cluster, and 
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Fig. 12. Ratio of the observed (extinction corrected) to the model Ha luminosity vs. the age of the clusters (left-panel), the mass 
of the clusters (central-panel), and the galactocentric distance for clusters with log(age) < 6.5 (right-panel). In these plots, the 
model Ha luminosities of clusters with masses M < 1450 M , were determined with the simulations described in the text, taking 
into account the incomplete sampling of the IMF. The error bars include both the dispersion of the model Ha luminosity and the 
uncertainty due to the error on the age of the clusters. Symbols are the same as for Fig. [4] 



it will be steeper. Assessing the details of the IMF in M33 and its 
variance requires the availability of deep colour-magnitude dia- 
grams for the stellar members of each cluster. This is difficult to 
attain in compact young clusters and is clearly beyond the scope 
of this study. 

8. The gaseous environments of young clusters 

Finally, we investigated whether the mass of young clusters is 
related to the surface density of the different phases of the sur- 
rounding ISM (ionised, atomic, molecular gas). 

As described in Sect. 2.4, the spatial resolution of the CO 
and Hi maps is 13"and 24"x 48 "respectively. Thus, the size of 
the clusters (see Table [2]) is smaller than the resolution of the 
CO images for roughly half of the objects in our sample, while 
none of the clusters can be resolved with the beam size of the 
21 -cm map. For this reason we limit our analysis to areas of the 
ISM that are larger than the clusters. Using an aperture equal to 
the lowest resolution gas map (Hi) i.e. of 24 arcsec, centred on 
the optical position of each object, we measured the gas surface 
densities associated with it. 

Figure [13] shows the measured surface density of the differ- 
ent gaseous components as well as the total gas surface density 
compared to the cluster stellar mass. There is no trend with clus- 
ter age, the surface density of the surrounding ISM is not affected 
by the early stages of cluster evolution. The ionised gas seems 
to better correlate with the cluster mass (as expected, since the 
most massive clusters have a higher number of ionising pho- 
tons). Concerning the Hi gas, one can only infer from the fig- 
ure that the majority of the clusters are found in a narrow re- 
gions of surface density (Nhi = 1 -r 3 x 10 21 cm -2 ), which is no 
big surprise given the rather flat Hi distribution in the M33 disc. 
Clusters are born in correspondence to Hi filaments, but there is 
no clear connection betwee n the stellar mass and the gas surface 
density dVerlev et al.ll2010h . 

The bottom-right panel of Fig. [13] shows the sum of all 
gaseous components (ionised, atomic, molecular) versus the 
stellar mass. It shows that the gaseous environment of these clus- 
ters is dominated by the atomic component. For a change of two 
orders of magnitude in cluster mass there is only half an order of 



magnitude change in the total gas column density. Because the 
scatter of this correlation is large, as Verley et al. (2010) have 
already pointed out, the gas surface density is not the only pa- 
rameter which determines the star-formation rate and the mass 
of the newly born clusters. 



5.0 

^ 4 - 5 

2 4.0 

i 3.5 

CP 

5 3.0 

2.5 
2.0 

5.0 

4.5 

<^ 4.0 

5- 3.5 

3 3.0 

2.5 
2.0 



. HII 

> : 

• • 

• 


. HI 

V : 

* > • 

• 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

■H 2 

A 

•*: • 

• 

• 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

■ H T0T 

• 



-0.5 0.0 0.5 1.0 0.5 1.0 1.5 2.0 
log 2 [M e pc" 2 ] log S [M Q pc" 2 ] 
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codes are the same as for Fig. [4] 



9. Summary and conclusions 

We analysed multi-wavelength observations of a sample of 
young star clusters and associations in M33 combining UV 
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(GALEX), optical (UBVRI), and MIR (Spitzer/lRAC-MlPS), 3- 
mm (BIMA-FCRAO) and 21 -cm (WSRT) observations. The UV 
and optical data sets were used to derive their age and mass, and 
the line-of-sight extinction by fitting the measured SEDs to a 
set of single-burst stellar population models. The model SEDs 
were reddened using various extinction laws (Milky Way, LMC, 
SMC), and the best-fit solutions were obtained with the LMC 
laws. 

We find ages between 2 and 15 Myr and stellar masses be- 
tween 3xl0 2 and 4 x 10 4 M . The extinction towards these 
stellar systems, Ay varies between 0.3 and 1 mag, indicating 
that they are not highly obscured even in an early evolutionary 
phase. As expected, older clusters are less extinct than younger 
ones. We also find a correlation between mass and radius with 
Moc Ri-5±o.2_ There is no clear trend between extinction and 
other parameters such as metallicity or total gas surface density. 
The fractional luminosities at 8 /im and 24 /im decrease with 
the galactocentric distance, suggesting a lower PAHs and small- 
grain abundance in the outer regions of the galaxy. A comparison 
between the model SEDs and the Spitzer/IRAC fluxes shows an 
excess emission at 3.6 pm and 4.5 yum, which may be due to a hot 
dust component, heated by the ionising radiation of the cluster 
massive stars. 

We find a discrepancy between the observed Ha luminos- 
ity and that expected from population synthesis models, espe- 
cially in young (< 3 Myr) and low-mass (M < 3 x 10 3 M ) 
clusters. The ratio of the observed-to-model Ha luminosity in- 
creases with cluster age and ranges between 20% and 90%. We 
showed that this difference can be reduced but not removed if 
one takes into account the uncertainties of the age from the SED- 
fitting procedure and the stochasticity of the IMF at the high 
mass-end. Similarly, the gap between the best-fit bolometric lu- 
minosity (derived assuming a fully sampled IMF) and the total 
observed luminosity (UV and IR) cannot be fully explained if 
only a stochastically sampled IMF is considered. We then dis- 
cussed other possible scenarios: 

1 . Absorption of Lyman continuum photons by dust in Hn re- 
gions cannot fully explain this discrepancy because the IR 
emissions inferred from the 24 /im luminosity is generally 
lower than expected from the derived extinction. This sug- 
gests that UV continuum photons are missing, but not be- 
cause of dust absorption only. 

2. The escape of UV ionising and non-ionising photons from 
the Hn region because of a clumpy and porous ISM is a likely 
scenario that would be supported by the large UV and Ha 
diffuse fractions observed in this galaxy. 

3. Finally we discussed the possibility of a delayed formation 
of the most massive stars in these clusters as a consequence 
of the increase in the observed-to-expected Ha luminosity 
ratio with age and mass. Because the oldest clusters in our 
sample are also the most massive ones, a stronger suppres- 
sion of massive star formation in low-mass systems than 
what is predicted by a stochastically sampled IMF would 
also alleviate the problem. 

Finally, we analysed the surface density of the gas (ionised, 
atomic, and molecular) in a region about 200 pc wide around the 
cluster. A marginal correlation only is found between the cluster 
mass and the total gas surface density, implying that on these 
large scales this is not the only parameter that determines the 
total mass of stars. 
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